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Abstract. We present a new method to compute the deflection of light rays in a perturbed 
FLRW geometry. We exploit the properties of the Geodesic Light Cone (GLC) gauge where 
null rays propagate at constant angular coordinates irrespectively of the given (inhomoge¬ 
neous and/or anisotropic) geometry. The gravitational deflection of null geodesics can then 
be obtained, in any other gauge, simply by expressing the angular coordinates of the given 
gauge in terms of the GLC angular coordinates. We apply this method to the standard 
Poisson gauge, including scalar perturbations, and give the full result for the deflection effect 
in terms of the direction of observation and observed redshift up to second order, and up to 
third order for the leading lensing terms. We also compare our results with those presently 
available in the literature and, in particular, we provide a new non trivial check of a previous 
result on the luminosity-redshift relation up to second order in cosmological perturbation 
theory. 
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1 Introduction 

It is well-known that the light-like signals emitted by the astrophysical sources propagate 
along the null geodesics of the cosmic geometry, and that the observed properties of this 
radiation (such as wavelength, polarization, propagation direction) are affected by the geo¬ 
metric properties of the large-scale space-time. For instance, if the geometry is expanding, 
then the frequency of the received radiation is shifted with respect to the emitted frequency 
according to the well known cosmological redshift effect. 

Here we are interested, in particular, in the shift of direction of the received radiation 
with respect to the angular direction of the source, referred to a system of polar coordinates 
centered at the position of a static observer. More precisely, if we denote with a = 1,2, 
the angular coordinates of the source in the given reference frame, and with a = 1,2, 

the polar angles controlling - in the same frame ~ the direction of the received radiation 
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at the observer position, we are interested in computing the general relation 0“ = ^“(^ 0)1 
determined by the given model of cosmological geometry. 

If the geometry is spatially homogeneous and isotropic the above relation, of course, is 
trivial: 0“ = However, if the cosmic geometry deviates, even perturbatively, from exact 
homogeneity and isotropy (for instance, because of macroscopic fluctuations of primordial 
inflationary origin, or because of the properties of the local matter distribution), then the 
light-like signals are geometrically deflected, and we have in general a non-trivial relation 
0“ = dictated by the angular profile of the distorted shape of the null light-cone 

hypersurface. 

The explicit form of such an angular relation depends not only on the given geometry 
but also on the chosen gauge. The aim of this paper is to present a new approach to 
the computation of the above angular relation in a generic (inhomogeneous/anisotropic) 
cosmological geometry, and apply this new method, in particular, to a perturbed FLRW 
metric which includes scalar perturbations up to third order parametrized in the usual Poisson 
gauge. 

Our method is based on the coordinate transformation relating the angular coordinates 
of the source 0“, expressed in the Poisson gauge (PG), to the corresponding source coor¬ 
dinates 9g expressed in the so-called geodesic light-cone (GLC) gauge [1], namely on the 
transformation 0“ = 6^(9^). Indeed, in the GLC gauge, null geodesics are characterized by 
constant values of the angular coordinates [1, 2]: it follows that 9^ = 9^, i.e. that the angular 
position of the source, 9^, always coincides in the GLC gauge with its “apparent” position 
determined by the local direction of the light ray received by the observer. In addition, the 
coordinate transformation between PG and GLC gauge can always be defined in such a way 
that the angular coordinates of the two frames coincide (to all orders) at the observer position 
[2, 3], namely at the origin of the PG system of polar coordinates: this implies 9^ = 9^. As 
a consequence, we have 

0“ = Gr(^) = GT{9l) = GT{9l) , (1.1) 

where GT denotes the above-mentionad gauge transformation, with the dependence on r not 
explicitly shown. Eq. (1.1) provides the sought relation for the geometric deflection effect 
induced by the cosmic gravitational field, parametrized in the PG. 

We recall that there is a large literature on the perturbative evaluation of cosmologi¬ 
cal observables, with and without the use of the GLC gauge. Using the GLC coordinates, 
expressions for the luminosity distance-redshift relation were obtained in [3-5] (with applica¬ 
tions discussed in [6-8] ^) and for the galaxy number counts in [10]. Using different methods, 
similar expressions were obtained in [11, 12] for the luminosity distance-redshift relation, and 
in [13-16] for the galaxy number counts. Furthermore, results for the Newtonian density 

^See [9] for further details about the application of the GLC coordinates in the framework of the light-cone 
averaging procedure. 


- 2 - 



fluctuation were given in [17] and the study of lensing up to second order was discussed 
in [18]. 

The paper is organized as follows. In Sect. 2 we briefly recall the main definitions and 
properties of the coordinate systems (GLC and PG) to be used in this paper. In Sect. 3 
we first evaluate the cosmological deflection effect starting from the transformation express¬ 
ing GLC metric and coordinates in terms of the PG ones and then inverting the obtained 
transformation, up to the desired order. We shall use, to this purpose, the perturbative 
second-order results already presented in [3], expressing them in terms of observational co¬ 
ordinates and extending them to third perturbative order (but only for the leading lensing 
contributions, i.e. for those with the maximum number of angular derivatives). In Sect. 4 
we present our new method based on the coordinate transformation directly expressing PG 
quantities in terms of GLC ones. Then, as an independent consistency check, we show that 
the results obtained for the deflection angle by this new method exactly coincide with those 
of Sect. 3. In Sect. 5 we show that our leading lensing terms satisfy, up to third order, 
a non-trivial “lens equation” and compare our results with those currently available in the 
literature through different computational methods. A few conclusive remarks are finally 
reported in Sect. 6. 

In Appendix A we provide another application (and an additional test) of our new 
method by computing, up to the second perturbative order and with the approach of Sect. 
4, the full result for the luminosity distance of the light source as a function of direction 
of observation and redshift, and show that it exactly coincides with the one obtained in 
[3, 4] using the approach of Sect. 3. In Appendix B we clarify the origin of an apparent 
disagreement on which averages of the luminosity distance are least affected by lensing. In 
Appendix C we prove that, in spite of its form, the result of Appendix A and of [3, 4] for the 
luminosity distance-redshift relation is actually covariant. 

2 Geodesic Light-Cone and Poisson gauges: a short reminder 

In this Section we shortly recall the main definitions and properties of the coordinate systems 
to be used in this paper. 

The so-called geodesic light-cone (GLC) coordinates [1] are particularly adapted to de¬ 
scribe signals that propagate along our past light-cone. They consist of a timelike coordinate 
r, of a null coordinate tc, and of two angular coordinates 0“ (a = 1, 2). The parameter r can 
be identified with synchronous gauge time [2], and thus provides the four-velocity of a static 
geodesic observer in the form = —d^r. The line-element of the GLC metric reads: 

ds^ = T^dw^-2TdwdT + -fab{dr-U'^dw)ide’’-U’^dw) , a,b = 1,2 , (2.1) 
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or, in matrix form: 


/ 0 -T 0 \ 


/ -1 

_y —1 
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( 2 . 2 ) 
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V-([/“)^/T 

0^ 
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where 0 = (0,0), Ub = (Ui, U 2 ) and U‘^ = jabU^'U'^. Here T, [/“ and jab = Jba are arbitrary 
functions of the GLC coordinates, and jab and its inverse 7 “^ lower and raise the two- 
dimensional indices. 

The w = constant hypersurfaces define a foliation of space-time in terms of the null 
hypersurfaces {d^wd^w = 0 ) that correspond to the past light-cones of a given observer 
throughout its history. Also, in this gauge, the null geodesics connecting sources and observer 
are characterized by the tangent vector = —ujg^'^du'w = = ujT~^Sr (where cu is an 

arbitrary normalization constant), meaning that photons travel at constant values of w and 
9°'. This crucial property of the GLC gauge will be extensively exploited in the following. It 
also renders the calculation of the redshift and luminosity distance (the two entries in the 
Hubble diagram) particularly simple. 

Denoting by subscripts “o” and “s”, respectively, quantities evaluated at the observer 
and source space-time position, one finds that the exact expression of the redshift Zs associ¬ 
ated with a light ray going from “s” to “o” is simply given (for static observers and sources) 


by [ 1 ] 

, ^ ^ ^ T{wo,To,6^) 

^ ^ {kf^U,)a {df^wd,r)a T{Wa,Ts,6-) ' ’ 

Namely, it factorizes in terms of an entry of the GLC metric evaluated at the observer and 
the same quantity evaluated at the source, precisely as in the case of the FLRW metric. 
Similarly, an exact factorized expression for the so-called Jacobi Map [19] was derived in [4]. 
The determinant of the associated Jacobi matrix allows to express the luminosity (and area) 
distance, di and dA, entirely in terms of jab and its derivatives at the observer, namely: 


d‘l = {1 + Zs)^d\ = 4(1 -h ZsY 




[det [ur^drj°‘^) 


(2.4) 


where 7 = det (jab)- The above result has been used to obtain non-perturbative expressions 
for weak leasing quantities such as magnification, convergence, shear and vorticity [ 20 ]. 

In order to connect our results to those usually presented in the literature we also conve¬ 
niently introduce the so-called Poisson gauge (PG) [21] (sometimes referred to, at first order, 
as the “Newtonian” or “longitudinal” gauge). Neglecting vector and tensor perturbations, 
the PG metric takes the form: 


ds"^ = a^{r]) — dg^ (1 -|- 2<1>) -|- (1 — 24') [dr^ + r^d9^ + sin^ 9d(j)^) 


(2.5) 


Here a{r]) is the scale factor, g is the conformal time, and the (generalized) Bardeen potentials 
$( 7 , r, 0 “) and '^{g,r,9°‘) describe the first-order {(f, ijj), second-order i/;G)) and third 
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order scalar perturbations of a conformally flat FLRW metric background: 

Cb = ^ + + ^^(3), ^ (2.6) 

2 6 2 b 

Assuming the absence of sources with anisotropic stresses we shall put everywhere (p = ip, so 
that, to first order, = 'k = •!/). 

For the application discussed in Sec. IV it is also useful to redefine the Poisson Gauge 
using the light-cone-like (but not exactly null) coordinate = r] + r, so that the PG metric 
(2.5) takes the alternative form: 


ds'^ = g^^dy^dy’^ = o?‘{g) — 2drf' (<I> -|- ^') -|- (1 — 2'!') _ 2dgdg^^ 

-|- (1 — 2'!') (ry’*' — ry)^ (dd^ + sin^ 9d4>‘^) 


(2.7) 


The inversion of Eq. 
order, then gives: 


(2.7), including the contribution of scalar perturbations up to third 


LIU 

dpG 


/-1 + A 
-1 + A 


-1 +A 
B 
0 ^ 



( 2 . 8 ) 


where 


A = 2ip - Aijj"^ + 

3 

B = Alp + ip^"^^ + — Aipip^'^'^ + 16ip^ + Aipip^‘^\ 

7pg = {l + 2ip + Aip'^ + ip^"^^ + 8ip^ + Aipip'^'^^ + 7o^ 


7 o^ = - d) diag ( 1 , sin ^ 0'j , 


(2.9) 


and where the scalar fluctuations ip,ip^‘^^ ,^p) ,■ 0 ( 3 )^ (^(3) functions of the PG coordinates 
ry, , 6°‘. 

We are now ready to carry out the connection between coordinates and metric in these 
different frames, and to compute the deflection angles in the PG. In the next two Sections we 
shall adopt two quite different procedures: the first is the one followed in our previous papers 
[1-5]; the second is conceptually (and also, to some extent, technically) simpler. In spite of 
being very different the two procedures will be shown to give exactly equivalent results. 


3 Expressing GLC quantities in terms of PG coordinates and physical 
observables 

In order to evaluate the light deflection induced by scalar cosmological perturbations we 
will consider, in this section, the transformation relating GLG coordinates and metric to 
those of the Poisson gauge, and then express the result in terms of the observer’s detection 
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angle and redshift. We shall use the approach already investigated to second order in [3, 4] 
(for vanishing anisotropic stress) and in [5] (for non-vanishing anisotropic stress), extending 
the previous computations in order to describe the light deflection in terms of observational 
variables and to include the leading lensing contributions up to third order. 


3.1 The coordinate transformation 

Let us consider the general (exact) transformation connecting the (inverse of the) GLC metric 
(2.2) to the (inverse of the) PG metric (2.5), 

= > (3-1) 

where we have denoted by = (r/, r, 0“) the PG coordinates and by = (r, w, 0°') the GLG 
ones. Solving this set of differential equations for the variables (with the boundary 

conditions that the transformation is non singular, and that 0“ = 0“ at the observer position 
r = 0,7? = r/o), we obtain [3, 4], up to second order in the scalar perturbations: 

T = (3.2) 


with 


with 


with 


= ( f dr]'a{r]') \ , = a{r])P{r],r,9°') 

yJvin / 


r( 2 ) = 


^(2) _ ^2 {drPf + -it daP d^P (r/', r, 0“) ; 


' Vin 


w = 


= 7/+, 

^( 2 ) = 1 r dx 


u;(i) = Q(r/+,r/_,6»“), 

^(2) ^(2) 4-ijjd^Q + Q^Q Q^Q 9 °'); 

0a _ 0a(O) _|_ 0a(l) ^a(2) 


9“(o) 

Qa(2) 


0 a(l) = 1 r dx 
'Ho 


itdbQ {v+,x,9'^) 


(3.3) 

(3.4) 


(3.5) 

(3.6) 


rri- 


dx 


'•no 


(7y+,x,r). 


(3.7) 


We have defined -ft = diag(r ^,r ^ sin ^0), and we recall that the lower integration limit 
din represents an early enough time when the perturbations (or better their integrands) 
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were negligible. We have also used the zeroth-order light-cone variables = r/ ± r, with 
corresponding partial derivatives: 


dr, = d+ + d- , dr = d+ — d- , d± = 


dri± ~ 2^^^^ 


(3.8) 


Finally, we have defined 

P{v,r,9°')=[ dr]'^^'^'il;{r]',r,9° 


Q{v+,V-,9‘")= dx'il^{r}+,x,9°-) . (3.9) 


'Vo 


Using the above results we can then compute the non-trivial entries of the (inverse) 
GLC metric and we find, up to second order [3, 4]: 


'Y'— 1 _ 


a{r]) _ 


1 -|- d+Q — drP + drjW^'^^ H - {d-q — dr)T^'^^ — + 'i/’^ 


-V' {d+Q + drP) - drPd+Q - Jq daPdbQ 

a a a 

+ {d+Q - drP) , 

j-b = „-2 l^ab ^ 2^P) + + (a 6)] + 7^"' + 4^’^) - dr^^^^^dr^d’’^^^ 

+drr^^^dr9^^^^ + 2^ U^^dc9^^^'^ + (a o 6)1 + -ifdc9’^^^^dd9'^^^'^ 


(3.10) 


(3.11) 


-h 


■^-40^2) + ^ 5)J I, 


(3.12) 


3.2 Light deflection to second order 


As already stressed in the Introduction, we need now to express the angular position of the 
source, given by the PG coordinates as a function of the corresponding GLG coordinates 
9g which exactly identify the “observed” position of the source, 0“, determined by the di¬ 
rection of the received radiation as measured by a static PG observer (see Eq. (1.1)). We 
also want to give the hnal result for 0“ in terms of observational variables only: the above 
mentioned angular direction and the measured redshift of the source, Zg. 

In order to complete the first step one has to invert the coordinate transformation given 
in Eq. (3.7). This was done, for the first time and up to second order, in [3]. The result we 
obtained. 


go. _ ga{0) ga{l) , ga{2) _ ga _ 

2 


Jrio Jrjo 


dxjS’' dy db'ip{r]+,y,9°-) 


-h- 


1 


rv- 


dx 7 o / dy dbipiy+,y,9'^) 


L-JVo 

■V- 


'Vo 


dr. 


■ pv- 

/ dx 7^ / 

-'Jrjo Jrjo 


dx-f^'^ dy dd'iljiy+,y,9‘^ 


dx [7o7c + V' r + A“] (7+, X, 0“) 


(3.13) 


'Vo 
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(see Eq. (3.30) of [3]), can be conveniently expressed in conformal time rj, and written for 
the angular coordinate of the source as 0“ = + 9s^‘^\ with = 0“ = and 

with 


rvo rvo 

= -2 / d7?'7o“' / dr^"d,i,{r ^"), 

Jrjs Jri' 


6/“(2) =43 


rVo rrjo 

Jn, Jn' 


rVo 


rVo rVo 

dv'lo^dd / dr]' 


J J7?3 


/? 7 ' 


where 


+2 / d,?' [7o“^Cc(r?') + V’(r?')e“(V) + A“(7')] , 

'Jr]s 


1 f r 

Cc(?/') = “4 y + V’o-V'-2 J dr]"'dr^n,'ijj{r]" 

+H'’da 


rVo \ f mo 

dri"i,{ri") )dJ dv'"Hri'") 


f JjO 


f 


V^Vo 
rVo 


rvo 


r(V) = -jim (70 (d)'9fe I dr]”'ipir]”) j - 7o (^O^b J ^ dr]”'ilj{r]") 

mo 


rvo 


dv"lodb 


in' 


V’o - 2 / dr] d^iiiilr{r] ) 


/ rj" 


A“(?/0 = -27o''5c dr]"ijj{r]")] j^ dr]"da 


+ 


rvo 


V’o - '0(7 ) - 2 / dr]"drjiiilr{r]") 
J T}' 


lo'db { / dr]"''ll;{r]"') 


rVo 


rVo 


7o“'' 5, / d7'V(7")- 

Jn' 


(3.14) 


(3.15) 


(3.16) 


(3.17) 


(3.18) 


Note that, in the above equations and hereafter, partial derivatives like da are made with 
respect to 9°“ = 9^. Also, for brevity, we only indicate the conformal time r] as argument on 
the integrands over the line-of-sight, suppressing other arguments like 

According to Eq. (3.14), we can also introduce a so-called (first order) leasing potential 
(see for instance [22]), defined by 


^{r]) 


2 

7o - 7 



dv'- -^'0(70 


7o - 7 



(3.19) 


and such that, to first order, 7“ = 'jQ^db^{r]s), where 7 q^ = r^ 7 g^. 

In view of possible phenomenological applications it may be appropriate, at this point, 
to evaluate the terms providing the leading leasing contributions at large enough redshift 
(z > 0.3), to second order in the deflection angle. Let us note, to this purpose, that the 
corrections arising from writing the deflection angles in terms of the observed redshift (see 
below) are subleading, hence they will be neglected. Considering Eqs. (3.15)-(3.18), let 
us hrst select those terms contributing to and containing the maximum number of 

transverse derivatives (namely, three angular derivatives). Denoting such terms with we 
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have: 


na(2) _ . 


rvo 


dvio^dd I dr]''Tp{r]") 


rvo 


'Vs 


dr 


rVo , rVo 

/ dr ]'/ dr]"'il){r]") 

'Vs Jv' 


7o''dbf [ d'n"'ilj{T]")\ dc( [ dri"-f^'^db [ dr]'"Tp{r]'‘ 

\Jv' / \Jv' Jv" 

+ \lo^db drj"dc (^J dTj"'ipir]"')^ ■ (3-20) 

After some algebraic manipulations and integrations by parts one finds that the above equa¬ 
tion can be rewritten in a somewhat more useful form as: 


dr]'- - ^ 7 o‘'dcdb'ip{r]') 


Vo - Vs 


'Vs 


Vo - V‘ 


+ 


Vo - Vs 
rVo 




+ 


JVs 

dv' 


Vo - V 


db^iv) 


>vs iVo-v') Jv' {Vo-V") 

^lo'{ddlt')dc(jdv^'^Piv"')] db(Jdr("i}{r(") 


+ {dbT^^)t^drS^{vs)ddJJ^ij]s) 


(3.21) 


Finally, neglecting the subleading terms with only two angular derivatives acting on the 
metric perturbations, one obtains for the second-order leading contributions the following 
simple expression: 

(6'“^^^) Leading =-^— [ dv'- - dd^ {v')dcdh'4’{v) 

Vo - Vs Jris Vo - V 

= -^ r dry'^ (3 22) 

Vo - Vs Jr,s Vo - V 

This is in perfect agreement with the recent result obtained in [23]. 

Let us now express the full second-order result for 0“ in terms of the redshift Zg ■ Follow¬ 
ing the procedure introduced in previous papers [2, 3, 5] we consider the convenient expansion 
Vs = Vs^^ + "Vs^^ + • • • and -|- + • • •, where the zeroth-order parameters satisfy 

the condition ryi*^^ -|- = Wo = Vo and are defined in such a way that the observed redshift 

Zg is given by 

1 + "- = ^- P.23) 

a{Vs ) 

By expanding around vi^\ the exact expression (2.3) for Zg, using Eq. (3.10) for 
and the definition (3.23) for Zg, and expanding to first order also the light-cone condition 
w = Wo = const, one finds that the first-order corrections vi^\ are given by [3, 5]: 


vi^^ = 


Jig 


rVo 


V\\o - V\\g + Ipo - Ips - 2 / dv'dr^'-ip {v') 


_( 0 ) 

Vs 


/ Vo 

dv'ip [v') , 


_{ 0 ) 

Vs 


(3.24) 
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where we have defined T-Lg = {d„a/a)_{o) and we recall that nii^ = {drP){r]i^\rf'\o), with P 

'Hs '' 

given by Eq. (3.9), is the first-order velocity perturbation projected along the radial direction 
connecting source and observer (see also [4]). 

The desired expression for 0“ in terms of the observation variables 0“ and Zg can now 
be obtained by Taylor expanding the coordinate transformation (3.13), up to second order, 
as follows: 


0“(z„0,) = + + 

+ #(a,0“)W(#,fW,0„) + f(i)(a,0“)W(#,fio),^^^ . (3.26) 

Here 0“1^) and 0“1^1 are given, respectively, by Eqs. (3.14) and (3.15) evaluated at 0“. 

The time parameter in turn, is directly related to the observed redshift by Eq. (3.23), 
namely by the standard relation dzg = —(1 + Zg)'Hgdf]'^\ 

The first order part of Eq. (3.26) then gives: 

rvo 

ef^\zg,eo) = -2 d'n'jS^ dr]” dbi^ivi") (3.27) 

Jr)s Jv' 

(see also Eq. (3.14)). Considering the second-order part of Eq. (3.26), and keeping all 
(leading as well as non leading) terms, we obtain: 


O':^^\zg,9o 


J Full 


= 4 db 


rvo 


, -( 0 ) 


rVo 


dvio^dc / df] rpir] ) 


rVo 


dv'lo^dd I dr]”rp{r]”) 


rvo 


, -( 0 ) 
J Jrjl 




Vo\\ - I'sll + V’o - V’s - 2 


rVo 

/ , 0 ) dr]'dr^frP{r]') 
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rvo 


/-(O) 


dr]''il>{r]'] 
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rvo 


iVo-Vs) Jvi' 


( 0 ) 


dr]'%%rl;{r]') 


+4 


rVo rVo 1 rrjo 

dr]'ir{r]') dr]'- -^ / dr]''^S’^dbip{v”) 

Jf]i^ Jnr iVo-V) Jri' 


rvo 


rvo 


-4 / dr]' rpir]') / dr]'-fodb 


-{ 0 ) 

'Vs 

rvo 


-( 0 ) 

Vs 


V'(V) - V’o , r° J // a ,r 
----h / dr] dr,’>rlr{r] ) 


+2 / dr]' [7o“Cc(0') + ^iv')eiv') + A“(7')] 




(3.28) 


where are given by Eqs.(3.16)-(3.18). Hereafter we shall omit the terms from the 

velocity and the gravitational potential at the observer position, i.e. U||o and i/'o; since, in 
any case, the first can be removed by going to the CMB frame and, in general, they cannot 
be evaluated within perturbation theory. 

After some manipulations and integration by parts we can finally rewrite our complete 
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second-order result in the following simpler form 
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dri''^Q^db / drj” < -\- — 4'0(r/^O 
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-2 
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Vs 

rVo 


rVo rVo rVo 

d'n'-lfdbipi'n') / dr]''ilj{r]'') - 2 / dr]'i^{r]')ltdb / d?/'V(V') 

j)7' jf?r' j??' 


dri'itdb 


'vf^ 


rVo 


4{r]') + 2 / dvi"dr^ii'il){r]'') / dr]”'tjj{r]”) 


rVo 


' T)' / Jv' 

where is given by Eq.(3.21) evaluated for rjs = and 


(3.29) 


3.3 Leading third-order contributions 

The above computation will now be extended to describe the angular deflection of light-like 
signals up to the third perturbative order. We are interested, in particular, in the leading 
lensing contribution corresponding, at third order, to angular corrections generated by five 
angular derivatives of the perturbed cosmological metric. 

It will be enough to that purpose to consider only the leading terms of the coordinate 
transformation introduced in Sect. 3.1 for w and Let us recall here, for convenience, that 
such terms are given by 


yjio) ^ ^ ^{1) ^ Ql^r]+,r]-,d°-), 


( 2 ) ^ 


1 n- 


dx 


'Vo 


lo'^daQdbQ {r]+,x,9° 


and by 


QaiO) ^ Qa Qail) T 

2 4 

1 n 


lo’^^bQ {v+,x,9‘^), 


= dx {ri+,y,91- 


(3.30) 

(3.31) 

(3.32) 


'Vo 


Extending to third order the above results we obtain the following leading angular contribu¬ 
tions: 


0a{3) = ^ r dx 


'Vo 


where 


1 rv- 

= - dx 

° J Vo 


^bc 

7o 


7 o^ daQ 


+ dbQ 4 ^( 2 )^ + ^ab g^^(3)l (3 33 ) 


{v+,x,9^)db j 

J Vo 


dy i4^dcQddQ {r]+,y,9^) 


(3.34) 
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As before, we need to invert the above transformation in order to obtain 0“ as a function 
of the observation angle 0“ = However, as already mentioned, there is no need of 
performing also the redshift expansion (as done in the previous section) since that expansion 
would only produce subleading contributions. 

By inverting and Taylor expanding Eq. (3.33), order by order, around 0“ we then obtain, 
for the angular coordinate of the source. 


(0“^"^)Leading = dbOf^^ - 


(3.35) 

where all terms on the right-hand side are expressed in terms of Inserting the above 
results (3.31)-(3.34), and performing a series of algebraic manipulations, we eventually arrive 
at the following simple expression: 


(0^(3)) 


Leading 


T] 


Vo 


2 

- / dv' 

- Vs Jris Vo 


Vs ^ ab 


(6'‘'^^^)Leading(V)'9c5feV’(V) 


+ {v')dcdddbip{v) 


(3.36) 


where the first and second-order leading expressions for 0 “( 2 ) given by Eqns. (3.14) 
and (3.22), respectively. 

Let us conclude this section by including the contribution of the Bardeen potentials to 
second and third order in the calculation of the leading lensing contribution to This can 
be simply done substituting -0 with if; + (1/4) (?/(^) -|- -|- (1/12) inside the 

leading result of Eqs. (3.14), (3.22) and (3.36), and expanding up to the third perturbative 
order. We then easily obtain the following full leading results at second and third order 




1 1 


Full Leading 


fV 


2r]o-Vs 

2 


rvo 

/ drj' 

I ns Vo - V 

Vo 






drj' 

Vo - Vs Jns Vo - V 


"^-Xlodd^^ivWb^iv') 


(3.37) 


and 




1 1 


rvo 


Full Leading 


dv- - ^iS^db ip^^\v') + 


QVo-Vs 
1 1 
2r]o-Vs 

Vo - Vs 


Jns Vo-V 

rVo 

/ dra¬ 
ins Vo - V 

•Vo 


Full Leading {v')dcdbil{v') 


. ,V - Vs sab 

dv Z -7770 


'Vs 


Vo - V 


+ ^0''^^Hv)d'^^^Kv')dcdddb'ilj{v') 


(3.38) 


Eurthermore, anisotropic stress can be taken into account in this limit by simply replacing 
'll; with (?/ -|- </>)/2 as discussed in [5]. 
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We conclude this Section by recalling that we have neglected tensor and vector per¬ 
turbations, in despite of the fact they are inevitably generated from scalar perturbations 
at second (and higher) order. The lensing associated with tensor and vector perturbations 
induced at second order was already discussed in [3, 24]. In our context we want to underline 
that their contribution to lensing is subleading with respect to the scalar one since it always 
involves, order by order, a smaller number of transverse derivatives. In other words, our final 
results for the leading lensing terms are unaffected by vector and tensor perturbations while 
the full results are not. 


4 Expressing PG quantities directly in terms of GLC coordinates: a more 
direct approach 

We will now present a different approach, based on the direct computation of the PG angles 
0“ (and possibly other PG quantities) as functions of the observation angles 0“ (which, as 
already stressed, coincide with the angular coordinates 0“ of the GLG frame) and of the r 
coordinate of the source. An important advantage of this new method is that, in order to 
give the result in terms of just physical observables, it is sufficient to express in the final 
expressions in terms of redshift and GLC angles. Given the novelty of this method, our 
discussion here will be more detailed than the one in Sect. 3. 


4.1 The coordinate transformation 


We start with the transformation expressing the PG metric (2.8) in terms of the GLG metric 
(2.2), and given in general by 



dy^ dy'^ 
dx^ dx^ 



(4.1) 


where we have denoted by y^ = (r/, r/"'', 0“) the PG coordinates and by the 

GLG ones. We also recall that, if we neglect higher order corrections, the standard FLRW 
geometry is represented in the GLG parametrization by the following zeroth-order metric 
components [2, 3]: 

(y-l)(0) ^ jja{0) ^ ^afe(O) ^ (4 ^) 

and that, to this order, the GLG and PG coordinates are connected according to Eq. (4.1) 
by the (almost trivial) transformation: 


^(0) = 




j^+(0) ^ Qa(0) ^ 


(4.3) 


Let us now include perturbations up to the desired order, and let us expand the PG 
coordinates y^ in terms of the GLC ones around the above lowest order expression as: 


y^. = (yM)(0)(^) + (yn^^Hx) + {yn^^\x) + • • • . 


(4.4) 




Consequently, we expand around all components of the PG metric (2.8) as 


dpciy) = + f {y^)^^\x) + • • • . (4.5) 

V J y=yW 

Finally, we perform similar expansions, up to the required order, for all components of the 
(inverse) GLG metric (2.2): 

T-^(a:) = (T-^)(°) + (T-i)(^) + • • • , 

?7“(x) = + + ••• , + ... (4.6) 

Here we are interested, in particnlar, in the transformation providing the angular co¬ 
ordinates of the PG frame computed up to second perturbative order. To this purpose we 
need the first-order transformation of all PG coordinates rj, rj~^ and Hence, let us start 
applying and expanding Eq. (4.1) by including, for the moment, only first-order corrections. 

In this approximation, the computation of the ^pQ component of Eq. (4.1) then gives 
the following differential eqnation for 

= —— — , (4.7) 

a a 

where we have defined H = drO-ja. In a similar way, the computation of the transformation 
relative to the ^pQ component of Eq. (4.1) gives the differential eqnation for 

(4.8) 

a 

and the transformation of the gpQ component of Eq. (4.1) gives the differential equation for 

0a{l). 

, (4.9) 

where we recall that db denotes partial derivative with respect to It should be noted that, 
according to the expansion (4.5), all components of the PG metric (such as a, "0) appearing 
here and in the subsequent equations are always functions of the unperturbed PG coordinates 
(y^)(*^), namely, considering the relation given in Eq. (4.3), they are functions of the GLG 
coordinates {t,w,9^). 

Gonsidering the same boundary conditions as those used in Sec. 3, integrating with 
respect to r the last two eqnations, and using the fact that w and 0“ are constant along a 
null geodesic, we easily obtain, along such a curve, 


7+(i) = 



dr' 

a(rO 




(4.10) 


and the first-order expression for the deflection angle in terms of GLG coordinates is then 
given by: 



dr” 

o(r") 


db'ip{T”,w,9). 


(4.11) 
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Let us also report, for later use, the first-order result for the conformal-time coordinate 
The general solution of Eq. (4.7) gives 


where 


|y(l) 





= w - -|- 


(4.12) 


(4.13) 


as can be easily checked by differentiating Eq. (4.12) and inserting the results into Eq. (4.7). 

In order to compute the full second-order angular transformation, , we need also the 
expressions for the first-order perturbations of the GLG metric (2.2). They can be obtained, 
as before, from the components of the general transformation (4.1). Gonsidering in particular 
the transformation of the gpQ component we obtain, to first order, 

(T-i)(i) = (4.14) 

a 

or, using Eq. (4.7): 


(T-i)(i) = - 

a a a 

Similarly, considering the transformation of the component, we obtain 


(4.15) 


or, using Eq. (4.9), 

[/a(l) ^ _ ^abQ^^+{l) _ 

Finally, from the g^Q component of Eq. (4.1): 



(4.16) 


(4.17) 


(4.18) 


We are now in the position of extending to second order the transformation of the 
angular coordinates. We must compute, first of all, the differential equation determining to 
second order the light-cone coordinate 77'*'. Considering the g'pQ component of Eq. (4.1), and 
including all second-order contributions, we obtain the differential condition: 


5 .^ 77 + 12 ) = ^^( 2 ) jp 2 H'ipg^^'> 

2a \ /2a a 

-2 77 ( 1 ) dri^ - -r7+(i)a^7/; - - , (4.19) 

a a 

where we have used Eq. (4.15) for (T“^)*^^\ and Eq. (4.8) for ■ By integrating 

with respect to r along a null geodesic, as before, and using the first-order results of Eqs. 
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(4.10)-(4.12) for r], rj^ and we thus obtain (hereafter, for simplicity, we omit to write the 
dependence of 70 and ^|J on the constant arguments w and 0 “): 



(4.20) 


Finally, considering the transformation of the g^Q component of Eq. (4.1), and using 
the results of Eqs. (4.15), (4.16) and (4.18) for the first-order components of the GLG metric, 
we obtain, to second-order: 




(4.21) 


By inserting the first-order results (4.10)-(4.12) for 7 , and the second-order result 
(4.20) for and integrating with respect to r, one immediately obtains from the above 

equation the second-order contribution to the sought transformation, expressing the PG 
angular coordinates 0“ in terms of the GLG coordinates: 

dr'{r',w,6’’). (4.22) 

4.2 Light deflection to second order 

In order to reproduce the results of the previous Section we need to express the PG coordi¬ 
nates of the source, in terms of both the observation angles 0 “ (which, as already stressed, 
exactly coincides with the GLG angles 0“) and of the observed redshift Zg- 

Let us start considering, as in Sect. 3.2, the leading second-order part of the coordinate 
transformation corresponding to the lensing contributions to generated by three angular 
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derivatives, and contained inside the first three terms of Eq. (4.21) Let us call them . 
Using Eqs. (4.10), (4.11) for and and considering the second and the last line of 
Eq. (4.20) for the leading terms of 77 "*“^^^ we obtain: 


0 


a(2) _ 
s,B ~ 



dr" 


a(r'") 

dr"" 

U, 1 ^ 

a{T""] 

) 

dr'" 

ddi 

a{T"') 

■° dr'" 

7), 


ddipir'") 


a{T"') 


(4.23) 


Let us now go from the GLC coordinates (r, w) to the zeroth-order PG coordinates 
(j^(o)^ ^+(0)) ijging Eq.(4.3). Being already at second order, we can then drop the suffix 
(0) and move to the standard PG coordinates ( 77 , r) considering the zeroth-order relation 
w = rj + r = Wo = rjo- 

After several algebraic manipulations and integrations by parts we can then write the 
above result in the equivalent, more useful form: 


^<^(2) _ _ 


+ 

+ 


Vo - Vs Jris 

2 r^° 


dv' Z - 7r7“<9<i ^{r]')dcdb'il;{r]') 


Vo - Vs Jns 

4 no 


dvj 


Vo - Vs Jns 

1 


dvj 


Vo - V 
f V' - Vs 
Vo - v' 
,v' - Vs 


{doTo)dd^^{v') - m^vl 


dc^ir]'] 


rvo 


+ 57 ;^ 


rvo 

+ 4: dr]' 

Jr}s 


-7o^^ { db I dv"ir{ri") 
v' 

rVo 


Vo - V 
rVo 

dc / dr]"'tp{r]") ] { db j drCirir]') 

J Tj' 


rvo 


dd / dr] rfrir] ) 
J rj' 


db [ dr]''ij:{r]'')\ (ddZo'') [ dr]”'yo''(dc [ dr]'"'4r{r]‘ 

Jn' 7^2 y„/ \ J„„ 


(4.24) 


where we have used Eq. (3.19) for the lensing potential. Neglecting all subleading terms 
without three angular derivatives acting on the metric fluctuations we are left only with the 
first-line contribution, and thus we recover for (ds^^^)Leading exactly the same result as the 
one given by the method of the previous section, see Eq. (3.22) (also in agreement with the 
results presented in [23]). 

^ As in Sec. 3, the contribution arising from writing the deflection angles in terms of the observed redshift 
are subleading with respect to the leading lensing contribution. 
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We also note that, by applying the following identity, 



(^dd J , 


(4.25) 


it can be explicitly checked that the result of Eq. (4.24) for 0“ exactly coincides with the 
result for given in Eq. (3.21), i.e. that 0“^^^ = 9^}"^. 

Let us now provide the full second-order result for 0“ not only in terms of the observation 
angle 0“ but also in terms of the observed redshift Zg, related to the coordinate Tg by the 
exact expression (2.3). To this purpose, similarly to what done in Sec. 3 and in previous 
papers [2, 3, 5], we expand Tg around the convenient values defined in such a way that 
the observed redshift Zg is given by^: 


T{to,Wo,9o) ^ a{To) 
T{Tg,Wo, 9 s) a(ri°^) 


(4.26) 


By setting Tg = + • • •, by using for T ^ the first-order result (4.15), and expanding 

to first order the above equation around fg^^ we easily obtain 


T(ro) 

T(r.) 


«(To) 


I + ^I - -k 




(4.27) 


from which, by applying Eq. (4.26), we have 


f(1) = - 

" H 




(4.28) 


where everything is now evaluated at Notice that, following Sect. 3.2, we have neglected 
already at this stage the perturbative contributions evaluated at the observer position. 

The final expression for 0“ in terms of the observation variables (9^,Zs) can now be 
obtained by Taylor expanding the angular coordinate transformation around fg^'^ . Including 
all contributions up to second order, using the fact that w and 0“ are constant along the null 
ray trajectory and that, to zeroth order, 0“^°^ = 0“ = 0“) we can write: 


9tizg,9o) = 0“ + 0“(i)(rW,u;,0„) + 0“(2)(Ti°),u;,0o) 
+ ri^)(5,0“)«(rio),u;,0o) + -- - . 


(4.29) 


Here 0“(^\ ^“(2) ^j-g given, respectively, by the expressions (4.11), (4.22) evaluated (like 
fi^^) at Tg = fg^\ The time parameter fs^\ in its turn, is directly related to the observed 

^Another computational advantage of the approach of this section is that we need to replace only one 
model-dependent parameter, Ts, with the observable quantity Zs, instead of the two parameters qs and 
appearing within the approach of Sect. 3. 
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redshift Zg and to the parameter of Sect. 3.2 by the standard relation /a = d'q'^^ = 
—dzs(l + Zs)~^T-Ls~^ (where l-Lg = a)_(o), as before). 

Vs 

By taking into account all contributions we can now perform some useful simplification 
in the second-order part of Eq. (4.29). Switching to the conformal time parameter 
defined above, and separating the first and second order contributions, we then find: 


rvo rvo 

z,) = -2 / / dri»dtHri”y, 

Jvs Jr}' 


(4.30) 


/nN Pr]o I 'fu I '/u I 'lu 

= efj _ 2 / dr?V(^') 7 o“' 9 fe / dr?'V(r?") - 2 / j ^^'^'(r/") 

2 


rvo 


rrio 


rvo 


-(0) 

'Vs 


-(0) 

'Vs 


oH. 

rvo 


rvo 


7o db / drj'tpir]') d, 


'Vs 


_( 0 ) 

rVs 


dm!a{r]')il} [r]',C{Vs,r]')) 


Vin 




/-(O) 

'Vs 

rVo 


j drj” -4: j dr]" 


rvo 


[d^ipir]")) / dr]"''4){r]" 


17}" 




lvi°^ 


dvio 


rVo 


2'ip{ri')db / dr]"ii{r]") - 


Tjo - r]' 


rVo rVo 

/ dr]"'ijj{ri")db / dr]"'ijj{rj") 
Iri' Jr]' 

/// / / ///> 


/(O) /, [dd'lS’’) j^ dr]''j^^dc j d?7"V(?7' 




ns 


tps-2 du 


'V 

rVo 


drj' 


'Vs 


rVo 


/ d77V(7') 

J Vs 


(4.31) 


We can immediately notice that the first order result exactly coincide with the one 
previously obtained in Sect. 3, see Eq. (3.27). In addition, using the fact that 
performing further simplifications by noticing, for instance, that 




-(0) 

rvs 


rvs Q ( 

dr]'a{rj')'il;{r]',({r]s,r]'))= drj'-^dr'ip {r]',rjs - t]o,9o) = v^g , (4.32) 

J 7}iji Jrjin 

dr]'dyjii{r]',w,6) = j (^-^ + dv'J^ = + j d'q'dn'ilr, (4.33) 


and integrating by part several times, it can be eventually shown that the result of Eq. 
(4.31) is exactly equivalent to the full expression for obtained with the complementary 
approach of Sect. 3, and already reported in Eq. (3.29). 


4.3 Leading third-order contributions 

Following the same procedure of Sect. 4.1 we now extend to third order the computation 
of the angular coordinate 0“ as a function of 6^, considering however only the contributions 
of the leading lensing terms with hve angular derivatives of the metric fluctuations. As 
already stressed, in this leading approximation there is no need of performing also the redshift 
expansion. Considering the component gpQ of the coordinate transformation (4.1) we then 
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obtain the following differential equation for 



Leading 




ab 


+a 

V / Leading 


7 


+ 7t(d^0a{2)\ 

d \ / Leading 


atm 9^^+w 

/ Leading 

_ (a6r?+(2)) 


Leading 


(4.34) 


We can use the results of the previous sections for 0“ and rj'^ to the first and second order, 
but we still have to compute and 

Let us start with 7 “^*-^^ By expressing the component g^Q of the transformation (4.1) 
up to second order, by using Eqs. (2.8) and (2.9), and keeping only leading lensing terms 
with four angular derivatives, we obtain: 



Leading 


V / Leading CL‘^ \ / Leading 

+ (4.35) 


The above equation (4.34) for (0“)(^1 can thus be rewritten as 



Leading 



.,ad 


Leading CL 


Leading 


+ (dce'^^^A 5fery+(i) 

ins V / Leading 


(4.36) 


/ o\ 

Let us finally compute (^?''')Lading- purpose we first need the leading, third-order 

contributions to the component g'pQ of the inverted PG metric (2.8). By using Eq. (2.9), and 
expanding to third order the perturbation -0 with respect to the PG coordinates (according 
to Eqs. (4.4), (4.5)), we obtain the leading terms 


{9 


PG / Leading 


4 



Leading 


4 dc^ + 

\ / Leading d^ 


(4.37) 


By computing the i^pQ component of the transformation (4.1), and using for ( 7 “^)^) the 
result (4.18), we then find that the differential equation for ry"*" with leading, third order 
contributions, is given by 


/ Leading d V / Leading 2d 

- dci!^ - 

2ci CL \ / Leading CL 

(4.38) 


By integrating with respect to r, and inserting the result into Eq. (4.36), we eventually arrive 
at the explicit form of the leading third-order contribution to the angle 0 “. 
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Switching to the conformal time coordinate we can then compare such a leading expres¬ 
sion with the analogous result obtained with the complementary approach of Sect. 3, and 
presented in Eq.(3.36). After several algebraic manipulations and integrations by parts we 
can show that the two expressions exactly coincide, thus confirming the equivalence and the 
validity of the two approaches illustrated in Sect. 3 and Sect. 4. 


5 Our leading-order results and and the so-called lens equation 


In this Section we would like to establish a connection between our (leading lensing) results 
and what is sometimes referred to as the (linearized) lens equation [25]. This will also 
facilitate the comparison of our results with those in the literature. Let us introduce an 
amplification matrix A'^ as the derivative of the angular coordinates of the source with 
respect to the angular direction of the light ray received at the observer’s position, namely: 


A 


a _ 

b — 



(5.1) 


Using this definition, the results of Sect. 3 and 4 provide all we need for evaluating 

We want to stress, however, that the above definition of is not fully satisfactory from 
the theoretical point of view. One obvious objection is that (5.1) depends on the coordinate 
system one is using. A better way would be to start from the Jacobi map connecting a 
suitably projected displacement vector at the source to the angle at the observer (see, for 
instance, [19, 20, 22, 26]) and then introduce the angular coordinate of the source by dividing 
the displacement vector by an unperturbed distance between source and observer. This 
definition has the advantage of being gauge independent but instead depends on the reference 
background model used to define We refer to the nice review paper [27] for a thorough 
discussion of all the ambiguities one encounters in defining such a matrix as well as the correct 
equation it is supposed to obey. 

Here, more modestly, we shall adopt the definition (5.1) and show that its leading lensing 
contributions satisfy, up to the third perturbative order, the above-mentioned (linearized) 
lens equation. Equivalently, by defining the deformation part of the amplification matrix 
by subtracting the zeroth order contribution, i.e. by introducing the convenient quantity 
^b~^b~ show that such a quantity satisfies the following equation [25]: 


n = 


Vo - Vs 
2 

Vo - Vs JTIs 




jvS^dM{v',Vo-v',0<^)At 


Vo - v 


dr ]'^^ 
Jns Vo - V 


%‘^dM{v',vo-v',0'^ 


Sh - 


(5.2) 


whose structure is clearly well suited for an iterative (or in some cases perhaps even an exact) 
solution. 
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In order to check that this lens equation reproduces our results up to third order we 
expand it up to the desired re-th order by setting = 0: 

Qa(n) 


(^a)(n) ^ _ 


dOs 


09^^ ’ 


n > 1. 


(5.3) 


We will obtain, in this way, iterative solutions for the deformation matrix and for 

the corresponding angular deflection Stopping the computation at first order we have 

[25], in particular, 

2 


(^»)(i) = 

which corresponds to 


dr]' -^ %^dcdbij{r]', r]o - rj', 9^), 

Vo - Vs Jris Vo - V 


= -- 


dv- - Vo-V ,02), 


(5.4) 


(5.5) 


Vo - Vs Jrts Vo - V 

and exactly coincides with the first-order result (see e.g. Eq. (3.27)). A second-order 
computation gives 


iV - Vs .ac 

, . - 1 7o 

Vo - Vs Jris Vo - V 


/ dr]‘ 


dcdbdd'ilr{v')9'^^^^ - 5c<9dV’(?/')^f 


which corresponds to 


W2) ^ __ 


Vo - Vs 


rVo 


'Vs 


dv- - 

Vo - V 


(5.6) 


(5.7) 


and exactly reproduces the second-order leading result of Eq. (3.22) (also in agreement with 
[23], as already stressed). Einally, at third order we have 




rvo 


Vo - Vs 


dr]' 


V 


'Vs 


Vo - V‘ 


Vs .ac 

J ^0 


dcdbdMr]')9'^^^^ + ^d,dbdMW)0^^^^9<^'> 




(5.8) 


which corresponds to 



2 


Vo - Vs 



dr] 


, v' -Vs 
Vo - v' 


A,ac 

7o 


dcdd^iv') 9'^^^^ + ^d,ddde^Pivy‘^^^'^9<^') , 


(5.9) 


and exactly coincides, once more, with the leading order result of Eq. (3.36). This may sug¬ 
gest that the leading lensing terms are correctly resumed by solving (5.2) non-perturbatively. 

Let us conclude this section with a comment. The lens equation as defined by Eq.(5.2) 
differs by a sign in the last term with respect to the equation used in [28]. As a consequence, 
the results given in [28] for to second and third order differ from the ones we obtained 
here. Furthermore, as a consequence of this different sign, the results obtained in [28] for the 
deformation part of the amplification matrix cannot be seen as derivatives of the deflection 
angles with respect to the observed ones. This is in contrast with our starting definition 
(5.1) and with what we have just obtained, but, given the above-mentioned ambiguities in 
the definition of the amplification matrix, this point needs further study. 
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6 Conclusions 


In this paper we have discussed various aspects of the propagation of light-like signals in 
perturbed cosmological backgrounds making ample use of the so-called geodesic light-cone 
(GLC) gauge first introduced in [1], 

In practice, we have focussed our attention on the deflection of light rays, showing that 
the effects of inhomogeneities can be described, in an arbitrary gauge and to the desired level 
of accuracy, by simply considering (to the same level of accuracy) the coordinate (gauge) 
transformation that connects the given gauge (in our case the Poisson Gauge (PG)) to those 
of the GLG, since in this latter the angular coordinates remain constant along null geodesics. 
We have also verified that, when the deflection is expressed in terms of physical quantities 
(such as the direction of observation and the observed redshift), exactly the same result is 
obtained whether one starts from the transformation leading from the PG to the GLG gauge 
or one proceeds in the opposite direction. 

The explicit calculation was carried out for the PG deflection angle in both methods up 
to the second perturbative order and up to third order for the leading lensing contribution. 
Given the involved structure of these calculations, finding exact agreement between the two 
procedures provides a very non trivial check of their correctedness. Our expressions for the 
PG deflection angle to second and third order ( given in Eqs. (3.22), (3.29) and (3.36)), 
represent one of the main results of this paper. In that spirit we have applied the same idea 
in order to double-check the calculations of the luminosity distance-redshift relation dL{9'^, z) 
(in the PG and to full second order) already present in the literature in [3, 4] and in [11, 12]. 
Such a check is important since the published results are not in full agreement with each other 
[5]. The result for (11 (as a function of observer angle and redshift) reported in Appendix A 
is in full agreement with the one obtained in [3, 4] using our previous approach 

Goncerning the leading-lensing expression for the amplification matrix in the PG at 
second and third perturbative order there appears to be (barring typos) some important 
discrepancies between our results reported in Sec. 5 and those reported in [28]. This could 
alter the conclusions reached in [28] about the higher order corrections that weak lensing 
could have on the GMB temperature and polarization anisotropies, but we prefer to leave 
this interesting issue to further study. 
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A Consistency check for cIa up to second order 


The luminosity-redshift relation dL{z,9g) has been derived in previous papers [3, 4] ® in the 
limit of vanishing anisotropic stress and up to second order in the Poisson gauge (see [5] for 
its generalization in the presence of anisotropic stress), starting from the GLG gauge and 
following the approach summarized in Sect. 3. On the other hand, other results about the 
luminosity-redshift relation have recently appeared [11, 12]. These do not seem to be in 
complete agreement (see [5]) with those of [3, 4]. We believe that arriving at a commonly 
accepted expression for the perturbed luminosity-redshift relation is of fundamental impor¬ 
tance in view of present (see, for example, [29, 30]) and future (see [31, 32]) high precision 
observations of cosmic large scale structure. 

To this aim we compute here the luminosity-redshift relation, in the limit of vanishing 
anisotropic stress and up to second order in the Poisson gauge, using the approach presented 
in Sect. 4, and show that the result exactly coincides with the one obtained in [3, 4] (for 
simplicity, we neglect gravitational potential and velocity terms at the observer position). 
This represents, in our opinion, a highly non-trivial check of the correctness of the result first 
obtained in [3, 4]. 

Starting from Eq. (2.4) we first notice that, if we neglect the gravitational potential 
and the peculiar velocity contributions at the observer position, the luminosity distance of a 
source with general redshift z is just given by (see also [1]): 


di 


(i + T 


(7)y- 

(sin 0)^/^ 


(A.l) 


Therefore, we need to evaluate the full expression of jab at the second perturbative order. 
This can be obtained easily by considering the metric transformation (4.1) written in the 
form: 


_ PC 

de<^ ' 

where daV^ are exactly the quantities evaluated in Sect. 4. Hence, the induced metric on 
the two-sphere (and consequently the angular/luminosity distance) is given immediately in 
terms of the GLG angles (i.e. the observed ones) thanks to this new approach. 

According to the procedure of Sect. 4 we need, to this purpose, the explicit form of 
the second order quantities ry+(2) as well as the corresponding first order 


®The results of [3, 4] for the luminosity distance-redshift relation are equal to each other, apart from 
gravitational potential and peculiar velocity terms calculated at the observer position. These terms are 
neglected here. 
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expressions. For and we can use the results already presented in Eqs. (4.20) and 
(4.22), while for the result is still to be computed. 

Considering the i^pQ component of the coordinate transformation (4.1), we find that 
7 ^( 2 ) niust satisfy the following differential equation: 

= - adrT]^'^^ - ^ + aH 

— da'(p9°'‘^^^ + aHtprj^^^ — — adri^r]^^^ , (A.3) 


where we have defined H = aja, and a dot denotes differentiation with respect to r. In 
analogy with the first order case, we then find the solution: 


7y(2) _ 


1 


a(r) 


dr' 




+dwri8) 

+8^9°'^^^ daV!^^^ + yiif 

—c?aV’ — 11:8^7]^^'^ — a8ripri^^'' fr', ^(r, t'), 9^ 


(A.4) 


which, after some algebraic manipulation, can be rewritten in a more useful form as: 


7^(2) = 


1 


a(T) 


dr' 




where the following relations have been used: 

8w r'), 0“) = a{T')^r]^^^\T',^{T, r'), 0“) + 2 V’('r', C(t, r'), 0“), 


(A.5) 


dr' 


a^r«(r',C(r,r'),r) = a(r')^ ^^(r', ^(r, r'), 0“) - ?(t, r ), 0“), 


dr 


- rjW (r>, C(r, r'), 0“) = -aHr,W(r', ^(r, r'), 0“) - V’(r', ^(r, r'), 0“) ■ 


(A.6) 


Similarly, the result for ry"*"*^^) given in Eq. (4.20) can also be written in the more compact 
form as 


jy + ( 2 ) _ 




a(r') [; 

_|-jy + (l)Q^7y+(l) _|_ 0“(1) — lip . 


(A.7) 
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Considering now the g^Q component of Eq. (4.1), and extending the previous first-order 
calculation of we can express in terms of known quantities as: 


7 


ab 


( 2 ) _ 7 “ 


ab 


^+( 1 ) _ 77^^) 

2 a dr'i/j dci^ -h -h - 4 V’ -^- 77 ;- 777 


77+(2) _ jn(2) 

_9 J_ L _L Q 

^ 7 +( 0 ) _ ^( 0 ) ^ 


V 


.+( 1 ) _ 71 ( 1 ) 


n 


'+( 0 ) _ 


( 1 ) \ 2 

J — A'lp aH + AaH — 
g\^) I g 


+(i) _ fjW 
+(o) _ 7 ( 0 ) 


-2aHg^^'> + a^ (h^-H 


7 ( 1 ) 


H 


- 2 — 7 ( 1 ) r(i)5c7;(^ 


^+( 0 )_^( 0 ) + 2 ^ ^ a 2 a 2 

,,ab ^,cb 


^,ab ^.cb ^,ac 


/Cf ^ " ' 


„ac^,bd 
a 2 


-^7§''5e0“(i) + ^9dr(i)a,r(i) + 2 -7g-7(i)5er(i) + 2 


fed 

IL 

i 2 


H 


be „-|-(l) _ „(1) 


o 

ad 


q2 ^+( 0 ) _ |y( 0 ) 


^!^Qd{l)d^Qa{l) _ “^^caQ^ghil) ^ :^Q^Qc{l) Q^Qb(l) ^ ']^Q^Qa{l) Q^Qb{l) 


a 7+luj _ fj(u) 


9d7o“ 




(A. 8 ) 


By applying Eq. (A.l) we eventually obtain the following second-order result for the lumi¬ 
nosity distance relation: 


= a 


^1/4 

£ 1 ^ 

-h 

+ 

II 

sini /2 71 

'- 7 ( 0 )) 

fi , , 7+(i)-7(- 

Y '^+7?+(0)-7« 


H + 2H‘^^ - 

7 ( 1 ) - 


^+^^ 7 ( 1 ) + JVer(l) 


jy+(l) _ 7^(1) 

7 +( 0 ) _ 7^(0) 


'z z 7 +(o) — 7 ( 0 ) 

^ Ve r(i) + I aH 7 ( 1 ) V. ^'^(i) + I V. r(i) + I V. r(2) 

2 ' 


+x (V 


' - \ 5rfr(i)ac0''(i) - ^ f 


Sin ( 


(A.9) 


(recall that Va denotes the covariant derivative on the two dimensional sphere with metric 
da^ = dO"^ + sin^ dd^"^). 
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The above result is written as a function of the GLC coordinates (r, w, 0“). For a clear 
comparison with the analogous result obtained in [3, 4] we need to express dL not only in 
terms of the observation angle 0“ = but also in terms of the redshift Zs, related to the 
coordinate r* of the given light source. To this purpose we will apply the procedure of 
Sect. 4, extending to second order the redshift expansion since, in the case of di, even the 
zeroth-order term turns out to be r-dependent, i.e. drd^^'^ / 0. 

By setting Ts = fs^'^ + + fP\ and by Taylor expanding Eq. (4.26) up to second 

order around fs^\ we obtain: 

^ 

a{n ) a{n ) 

(A.IO) 

s 

In order for Eq. (A. 10) to be valid, order by order, has to satisfy Eq. (4.28) while, at 

-( 2 ) 

second order, we obtain the following equation for r* \ 


f (2) = jL 
" H 


2 




H - 


f{i) 

' S 


-n(2) 


+ a(T-i) 


(A.ll) 


where the right-hand side has to be evaluated at r = Ts^\ Using 

adr (T“^) + dri) — -, (A.12) 

a a 

we obtain 


r(2) = _ 

" H 




- aHrPfi^ - 2 (ff - + a (r-‘) 


n(2) 


(A.13) 


( 2 ) 

To fully evaluate Tg ' from the above equation we thus need the second order expression 
of T“^. Proceeding as in Sect. 4, and considering the coordinates transformation for gpQ, 
we first obtain: 


^ a 

a a 


T(i) 


8.,") + ^_ 


^ - 28 ri/')? 8 : — - H 

a a a 


. , drtt n't drtt f^\ (drCL)'^ 

+4 V’— v^ ’ - 2 — r]^^> + 2 ^ ^ 

a a a 


WY - 


2 dh 




(A.14) 


By applying the first-order results for 17^)“ and (T i)^^^ given in Eqs. (4.15) and (4.17), 
the results for given in Eqs. (4.7), (4.8) and (4.9), and by using 
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Eq.(A.3), we are then lead to the following result: 


(T - “iir - ^dy, 

+ -dr,ipr]^^^ + -dyjipr]^^^'^ + ^ (dy^A^A 

a a a a \ / Za \ / 

+ ^ “ ^dyjA^^dy;r]^^^'> 

+ da - A^A + ^OaA^^dbA"^ + ^5a7?+(i)afe7?+(i) 

a \ /2a 2a 

— —daA^^db'q'^^^^ . (A.15) 

a 

The sought expression for di in terms of the observation variables (0“, Zs) can now be obtained 
by Taylor expanding Eq. (A.9) around fA'^ ■ Including all contributions up to second order, 
we can write 


dL{zs,w, do) = d^^ArAKw, ei) + d^l\TA\w, 0“) + rA^'dr dA\TA\w, B'^A 

+Al\fA\wX) + A^'^drdf{fA\wX) + \ [tAA^ dUf{fA\wXo) 


+TA'^drd^AAA\w,eA, 


(A.16) 


from which, after a long but straightforward calculation, we can write our final result for 
dL{zs,w,9A, as obtained in the new approach of Sect. 4, as: 

i+(i) 


dL{zs,w,9] 


df\zs,w,9^ 


= 1 - V’ + 


7 ]' 


fj+{0) _ yy(0) 


+ +E ^p- dy,X^^^ + dy,A^A 


rjn 


JJ\ 


+ (0) _ ^y(O) 


+ ^ (^H + ZH^ ~ ~ 'Aj 

-adXv^"^ -dy,7PX^^^ -da^p9<^^^^ + 

--TpVc 9^^^'^ + -aH A^^ Vc 9<^'> + - j ~ j Vc 9<^'> + - Vc 

2 2 2 ?7+(o) - r/(o) 2 

^V’ — dyjX^^^ + dwA^^ — aHA^A 


+ 


+ 


2aH^ (r/+(0) - r?(o)) 

j7+(l) _ 77 ( 1 ) 

-2'>P + ^, - X + 


H 


rj+m _ 7 ^( 0 ) h2 ( 7 ^+( 0 ) _ 7 ^( 0 )) ' 

+aH (“ + 1) A^^ + ^ - 2dyuA^^ + 
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.,cd 
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'll) — + duiV^^^ — aHr]^^^ 

H - ^a^drli - aH'il)'q^^^ - ipd^, 

+(9^V'r/(^) + d-uji/) + da'ilj 9 °‘^^'> - d^ ^ 

+ (dw'n~^^^^^ + aHrj^^^dw - 2 

da (7/+(^) - r?(i)) + ^dar]^^^dbV^^^ + ^aa7?+(i)9b7/+(i) 

-iS'daV^^^dbTJ+^^A , 


(A.17) 


where we have defined H = 1 — [ 0 / 7 ( 7 ;^^*^^ — 7/*^*^^)]“^ and, for simplicity, we have omitted the 
suffix s. 

For an explicit comparison of this result with the one obtained in [3, 4], we have to 
rewrite Eq. (A.17) using the conformal time 77 and the radial coordinate r (or, equivalently, 
the zero-order light-cone coordinates 'r]~^,r]~) used in [3, 4], Recalling that dfs^^/a = 

(see Sect. 4.2), we have 


drA{TP,'w,9o) = {drj - dr)A{'r]^^\ ri°\ 9^) , (A.18) 

d^aA{f^^\'w,9o) = drA{'qf\rf\9^a) > (A.19) 

for any given quantity A = A{f^\'w,9o). On the other hand, considering Eqs. (3.8) and 
(3.9), we also have 


,^+(1) = _Q ^ (A.20) 

r^d) = _p ^ (A. 21 ) 

=71)- d+Q , (A. 22 ) 

d-wV^^^ = —drP- (A.23) 


Eurthermore, we note that the result of [3, 4] are expressed in terms of the variables 9Ad 
and as given in Eq. (3.7). Therefore, we have to express OAd and in terms of 
such variables. Using Eq. (3.13), and equating Eqs. (3.26) and (4.29), we then obtain 
0 “(i) = —qAA for the first-order quantities, while, after some algebraic manipulation, we are 
lead to the following relation for the second order variables: 


Var( 2 ) 


_V„r(2) + ( 1 ) 0 ^ 06 ( 1 ) + + 9A^^Va 

+daQdJ'^^^^ + Qd+ - ^Va (QTo'^bQ) + Va {Plo’^dbQ) • (A.24) 


Inserting now Eqs. (A.5) and (A.7) into Eq. (A.17), and taking into account the results of 
Eqs. (A.18)-(A.24), after a long but straightforward calculation, we obtain that the lumi¬ 
nosity distance-redshift relation obtained with this new method exactly coincides with what 
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obtained in [3, 4], both at first and second order. For example, considering only the first 
order result and using Atj = rjo — fis^\ we obtain 


d^l\zs,w,e^) 

S^\zs,w,e^) 


-i^s + 


+(i) 

Vs 


+(0)“ " (0) 
Vs -Vs 




+ i V, of) + S 
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- 1 - 
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PsAv 
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V's + 2 / dv'dr^^ilj {v) 

Jrjl 


rVo 


+ 1^ / dvi^iv) - 


1 rvo _ „(o) 

-*■ / J / '/ \ ! I ( ! 


^V JA°'> 


dv 


^V Vo - V‘ 


-A 2 V’ iv 


(A.25) 


in full agreement not only with [3, 4] but also with [33, 34], 

We believe, in conclusion, that the above extremely non-trivial check fully confirms the 
correctness and solidity of our old result (obtained for the first time in [3]) for the luminosity 
distance-redshift relation up to second order in the Poisson gauge. 


B Leading lensing contribution to luminosity related observables 


There has been some debate in the literature about which luminosity related observables are 
free from the leading lensing corrections. In this Appendix we try to clarify this issue. 

We start by summarizing the results obtained in [3, 4] for the luminosity distance- 
redshift relation ^£,(0“, z^), up to second order in scalar perturbations and in the Poisson 
gauge. We will assume no anisotropic stress® and keep only the leading lensing terms with 
four angular derivatives. Starting from the result of [3], and writing everything explicitly as 
done in [5], we have, up to second order: 
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(B.l) 


(B.2) 


®As seen, anisotropic stress can be taken into account in this limit by simply replacing tj) with {ip + (j))/2 


as discussed in [5]. 
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where Arj = rjo — , and where A 2 is the two-dimensional angular Laplacian referred 

to the PG angles at the observer position 0“ = {9o,4>o), i-e. = A 2 = -|- cot 6odg^ + 

(sin 9o)~‘^d‘^ 




From (B.l) and (B.2) we can then easily obtain the general solution for any power of 


dL- 


[dl) (1) 

(djyw 
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( 1 ) 


= n- 


( 0 ) 


[dl) (2) 


^ r 1 A 
= n— + -{n-l) 
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( 1 )' 


( 0 ) 


( 0 ) 


^LJ- dr ^ ' v<v 

With a little bit of work the leading second order results can be re-written as a term propor¬ 
tional to the square of first order lensing plus total derivatives terms. More specifically: 


(dl 
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rVo 


dv"il;{v'') db 


rVo 


dv"ip{v") 


,(B.4) 


where ^{v) is the lensing potential defined in Eq.(3.19). 

We note that for n = 2 the lensing-square term drops out and we remain with just 
total derivative terms (see, also, [23, 35]). As a consequence, if we integrate the inverse flux 
^di/d^^ (also called “reciprocal magnification” in [35]) around the background sphere we 
obtain no leading lensing contribution to both first and second order [35]. 

The above results for dL were already obtained in [7] but the claim made there was that 
the flux <1> ~ d^'^, rather than its inverse, has no leading-lensing bias (see also [6]). This is 
in agreement with the argument given for the first time by Weinberg in [36] (see also [37]) 
that the average of the flux over the celestial sphere has no leading lensing contribution in a 
universe filled with stochastically isotropic and homogeneous perturbations. 

The resolution of this apparent conflict is actually quite simple. In fact, in [7] the 
definition of angular average includes a proper-area weight on the constant redshift sphere 
T,{wo,Zs)- As a result (see Eq. (2.17) in [7]): 


(d„^)('Wo, 2s) ^ {1 + 2s) 


JdS 


Q 

~dS~ 


JdS 


— (1 + Zs) 


—4 J 

7^ 


— (1 + Zg) 


-4 


d-TT 


A{Wo,Zs) 


(B.5) 


with 

A{wo,Zs)= [ d‘^9°'^ (B.6) 

JT,{Wo,Zs) 

where (...) is the light-cone average defined in [1] . A{wo, Zg) is the proper area of T,{wo, Zg) 
computed with its induced metric 7 ab, and expressed in terms of the coordinates (tCo, Zg) that 

^This light-cone average was introduced by extending to null hypersurfaces the gauge invariant average 
procedure for space-like domains defined in [38-40]. 
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identify the particular deformed 2-sphere Ti{wo, Zg) on which supernovae at fixed z lie. Since 
we have y/j ~ [4], saying that the average of the over the directions of observation 

has no leading lensing contribution is equivalent to saying that there is no-leading lensing 
contribution for the area-weighted, angle-averaged flux over the past 2-sphere of constant 
redshift 

The above result also shows that the stochastic average introduced in Eq. (6) of [35] is 
equivalent, at the leading lensing level, to the light-cone average used in [7]. This is because 
the lensing correction appearing in Eq. (6) of [35] is equivalent to the lensing correction acting 
on the induced metric of the deformed 2-sphere T,{wo, Zg). 

In conclusion, the observable that receives the smallest bias from lensing depends on 
the precise definition of averaging: we have taken the point of view that the number of 
sources in a given solid angle is proportional to the proper area subtended by that angle, but 
ultimately it’s physics - and the way the observations are made - that should decide which 
is the mathematical definition that should be compared with the data [37]. 


C Covariance of our expression for di 


The luminosity distance di has been computed up to second order (and in the Poisson gauge) 
both without [3, 4] and with [5] anisotropic stress. However, the expression for di obtained 
in [3, 4] contains terms which, naively, do not appear to be covariant under spatial rotations, 
a property required for a physical observable like di. 

Considering, for instance, Eqs. (3.4) and (3.5) of [3] (or the final result given in Ap¬ 
pendix A), we find that the suspicious terms appear in the following combination: 

- , (C.l) 

2 4 4 sm 0 V / 

with and given in Eq. (3.7). Hence, the question to address is whether the 

expression in Eq.(C.l) is covariant under rotations. Let us first note that, for such second- 
order terms, the angles 0“ can be safely identified with the observation angles 0“ without loss 
of generality. 

In order to prove the covariance of Eq.(C.l) let us assume that the combination 


0a(2) _ 0a{2) 


1 

2 


0d{i)d^0^a)^ 


(C.2) 


behaves like a true vector. In that case we have that 


2 2 sm 0 


- ^0~"(')5d (V,0“(i)) , (C.3) 


and it follows that the sum of the terms of Eq.(C.l) is given by —(l/2)Va0“(^) plus a term 
which is explicitly covariant (because 0“(^) is covariant). Hence, to prove that Eq.(C.l) is 
covariant we only have to show that our assumption that 0“(^) is a vector is true. 

®See also [41] for similar considerations. 
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To this purpose, using the results (3.7) and (3.13)-(3.19), we can write the “non trivially 
covariant” (NTC) part of as follows: 


ga{2) 


NTC = ) - 2 / 


rvo 


'Vs 


CVo 


dv”i>{v”) dd 7o , 


(C.4) 


which, after an integration by parts, becomes: 
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(C.5) 

To check the covariance of the above expression we can explicitly introduce covariant deriva¬ 
tives and rewrite the above equation as 
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where T^^ are the usual Christoffel symbols. On the other hand, by integrating by parts, it 
can be shown that 


rvo 




'Vs 


dVV(7?")) 4^(7') rS/ = T<^df • 


(C.7) 


This implies that the second line of Eq.(C. 6 ) is vanishing confirming that 0 “O) ig indeed a 
vector. 
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